ADIABATIC TRANSPORT OF BOSE-EINSTEIN CONDENSATE IN DOUBLE- 

AND TRIPLE- WELL TRAPS 
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By using a close similarity between multi-photon and tunneling population transfer schemes, we 
propose robust adiabatic methods for the transport of Bose-Einstein condensate (BEG) in double- 
and triple- well traps. The calculations within the mean-field approximation (Gross-Pitaevskii equa- 
tion) show that irreversible and complete transport takes place even in the presence of the non- 
linear effects caused by interaction between BEG atoms. The transfer is driven by adiabatic time- 
dependent monitoring the barriers and well depths. The proposed methods are universal and can 
be applied to a variety of systems and scenarios. 



The trapped Bose-Einstein condensate (BEG) is now 
widely recognized as a source of new fascinating physics 
and interesting cross-over with other areas [l], In par- 
ticular, a large attention is paid to dynamics of weakly 
bound condensates. Here two main systems are usually 
considered: multi-level BEG in a single potential well 
(see |3j] and references therein) and a single BEG in a 
multi-well trap 0, S, H, 0, i, Hi • In the former case, BEG 
components are formed by atoms in different hyperfine 
levels coupled by the resonant laser irradiation. In the 
second case, all BEG atoms are in the same state but 
the trap is separated by barriers into a set of potential 
wells and BEG dynamics reads as tunneling between the 
wells. This case is realized in multi-well traps Q, arrays 
of selectively addressable traps fio'l and optical lattices 

iEI. 

The multi-component and multi-well systems are quite 
similar and can be treated on the same footing as a BEG 
with (weakly) coupled fractions. In the former, the frac- 
tions coincide with BEG components while in the later 
with populations of the wells. Even for the constant cou- 
pling, dynamics of multi-fractional BEG is very reach and 
varies from different kinds of Josephson oscillations (in 
analogy with internal and external Josephson effects in 
weakly bound superconductors, see discussion in [Hj) to 
Mott insulator transitions and other transfer regimes per- 
tinent to electrons in a crystal lattice ll| . The dynamics 
strictly depends on the interplay between the interaction 
of BEG atoms and the coupling 0, 

For a constant coupling, BEG dynamics is mainly re- 
duced to oscillating fluxes of atoms between the frac- 
tions. That was a subject of intense investigation during 
the last decade 0, H B IS 0, S, H ■ At the same time, 
still a little was done for exploration of irreversible trans- 
port of BEG between its fractions when their populations 
do not oscillate but evolve by irreversible and controlled 
way. Specifically, BEG transport means that condensate, 
being initially in one of the fractions, is then completely 
transferred to another (target) fraction and kept there. 
Being realized, BEG transport could open intriguing per- 



spectives for investigation of topological states [ij, , 
generation and exploration of geometric phases (which 
may serve as promising information carriers [l3 | ) , various 
implementations in quantum computing (e.g. to build 
analogs of Josephson qubits in superconductors [l3|), etc. 

Nowadays adiabatic population transfer methods seem 
to be most promising to carry out the transport of 
BEG. A variety of suitable adiabatic methods was al- 
ready developed for the aims of atomic and molecular 
spectroscopy [l^, between them Stark Ghirped Rapid 
Adiabatic Passage (SGRAP) [itI and Stimulated Rapid 
Adiabatic Passage (STIRAP) [l8l|. These methods ex- 
ploit two-photon population transfer schemes and intend 
to excite non-dipole states in electronic systems. How- 
ever, the methods are quite general and can be upgraded 
to other systems and scenarios, see e.g. 

The aim of the present paper is to adapt the SGRAP 
and STIRAP ideas for BEG transport and develop the 
relevant population transfer schemes. For this aim the 
proper time-dependent protocols of the system parame- 
ters (well depths and barrier penetrabilities) will be pro- 
posed. We will concentrate on double- and triple-well 
traps where the necessary control of the parameters can 
be easily realized, e.g. the barriers can be monitored by 
varying depths and separations of the wells [^, ■ 

By our knowledge, the SGRAP has not still employed 
in BEG transport. Quite recently an alternative tech- 
nique, Rabi switch, was suggested, where time-dependent 
Rabi coupling of a fixed duration was used to transfer 
solitons and vortices Here we propose in principle 
another scheme where both Rabi frequency and detuning 
are adiabatically controlled. Such SGRAP-like scheme 
has the advantage to be rather insensitive to the process 
parameters, e.g. the coupling duration 1^ 12 1- 

As for STIRAP, it was already proposed to transfer 
individual atoms 21[ and BEG In BEG the nonlin- 



earity due to the interatomic interaction U was shown 
to be detrimental for the adiabatic transfer (like dynam- 
ical Stark shifts in multi-level system) 0]. As a remedy, 
a large detuning A of the well depths (trap asymme- 



try) was proposed. We will show that, unlike conclusio 
[7|, both large non-linearity and asymmetry are det 
mental for STIRAP but the transfer is robust at th( 
modest values. Instead, in SCRAP-like scheme a lar 
time-dependent detuning A{t) is useful and even cruci. 

Altogether, we will develop and justify simple and i 
fective adiabatic schemes for BEC transport in doubL 
triple- and multi-well traps. 

Our calculation have been performed in the mean-fie 
approximation by using the non-linear Schrodinger, 
Gross-Pitaevskii equation for BEC: 



thi>ir,t) = [-—V^+Ve.tir,t)+go\^ir,t)\^]-^if,t) ( 

where the dot means time derivative, ^'(r, t) is the ck 
sical order parameter of the system, Vext {r, t) is t 
external trap potential involving both (generally tiir 
dependent) confinement and coupling, go = 4Tra/m 
the parameter of interaction between BEC atoms, a 
the scattering length and m is the atomic mass. 

For BEC with M fractions, the order parameter can 
expanded as [1] 



M 



^if,t) = %/7V^^fe(t)$fc(f) 



(2) 



fe=i 



where $fe(r) is the static ground state solution of ([T]) 
for the isolated (without coupling) k-th well [23|] and 
ipk{t) = yiVfc(i)e*'^''(*) is the amplitude related with the 
relative population of the k-th well Nk (t) and correspond- 
ing phase 4>kit). The total number of atoms N is fixed: 

jdf\^{r,t)\yN = j:tLim)^i- ^ 

Being mainly interested in evolution of populations 
Nk (t) , we dispose by integration of the spatial distribu- 
tions $fc(r) and finally get 0, Q 



M 



where 

is the coupling between BEC fractions, 

hEkit)^ J dri^lVm^ + nVextrnk] 
is the potential depth, and 

hU^go f df\^k\^ 



(3) 



(4) 



(5) 



(6) 



labels the interaction between BEC atoms. The values 
flkj{t), Ek(t), and U have frequency dimension. 




FIG. 1: Schemes of the adiabatic population transfer in 
double-level a)-b) and double-well c)-d) systems, a) Double- 
level scheme with ofF-resonant pump laser pulse flp{t) and 
detuning A(t). b) Time evolution of pump Qp(t) and Stark 
i}st{t) laser pulses, c) Double- well scheme with barrier cou- 
pling i}{t) and detuning A(t). The condensate (dot at the 
left well) should be transferred to the right well, as indicated 
by the long arrow, d) Time evolution of the coupling and 
detuning. 



We use the coupling of the Gauss form with a common 
amplitude K: 

rikjit) = KQkjit), Ofc^(t) =exp{-^^M_^} (7) 

where t^j and T^j are centroid and width parameters. 
Then, dividing ^ by 1/2K, one gets 

1 

thiPk = [Ek{t) + A|V'fe|']i^fe - ^ ^ ^kAt)^j (8) 



where 



Ek{t)^Ek{t)/2K, A^UN/2K 



(9) 



and the time is scaled as 2Kt — ^ t and so is dimensionless. 
In IH]) the key parameter A is responsible for the inter- 
play between the coupling and interaction. In principle, 
one may further upgrade ([8]) by canonical transforma- 
tion of Nk and 4)k to other unknowns, population im- 
balancies and phase differences, and thus removing from 



dS]) the integral of motion A'' [2j| . It should be empha- 
sized that, unlike the previous studies of the oscillating 
BEC fluxes in traps with constant parameters 0, 0, S 0] , 
we will deal with irreversible BEC transport by monitor- 
ing time-dependent panameteis, namely depths Ek{t) and 
couplings flkjit). 

In the present study, the cases with M=2 and 3 are ex- 
plored. For M=2 the SCRAP scheme for a double-level 
system [lil . 17 [ and its counterpart proposed for BEC 



transport in a double-well system are illustrated in Fig. 



A=0 

A(t)=0.003(t-47) 



A(t)=0.3(t-47) 




FIG. 2: Populations iVi(t) (solid curve) and N2{t) (dotted 
curve) calculated with different interaction A and detuning 
A{t). Left panels: no interaction (A = 0) and different rates 
of detuning A{t). Right panels: the fixed rate of detuning, 
A(t) = 0.3(t — 47), and different interaction. The coupling 
f2i2(t) is characterized by ti2=47 and ri2=7.1. Time is di- 
mensionless. The initial phase difference is zero. 



1, see plots a)-b) and c)-d), respectively. In both cases, 
we deal with time-dependent detuning and coupling. In 
the double-level system, the detuning A{t) {— dynamical 
Stark shift) and the coupling are generated by the Stark 
flst{t) and pump Qp(t) pulses, respectively. Instead, in 
the double- well system the detuning A{t) = E2[t) — Ei (t) 
is produced by variation of well depths Ek (t) , while the 
coupling rt{t) is determined by tunneling ([4]) between the 
barriers. The later can be easily monitored by variation 
of the well separation. Obviously, the cases a) and c) 
are very similar: they involve two states |1) and |2) to- 
gether with the time-dependent coupling and detuning. 
And in both cases the adiabatic population transfer takes 
place when the detuning and coupling are simultaneously 
and slowly varied. To make the transfer irreversible, the 
system has to cross only once the resonance (symmetry) 
point. For this reason, the pump pulse should cover only 
one of the flanks (left in the plot b)) of the Stark pulse. 
The same aim is achieved in the double-well system by 
crossing once the symmetry point A{t) = during the 
detuning rise, see plot d). In both two-level and two- well 
cases the adiabatic condition lld| 



\9it)\«JnUt) + AHt) 



(10) 



has to be fulfilled, where 6 is the mixing angle of the 
states |1) and |2) during the adiabatic evolution. 

In Fig. 2 the results of our calculations for the double- 
well case are demonstrated. The plots a)-c) show that 
the population transfer improves with increasing the de- 
tuning rate A(t) from 0.003 to 0.3. For A(t)=0.3 the 
transfer is complete. So, the robust transport of BEG 
needs a rapid change of the detuning, which agrees with 
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d) 

(2) 
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FIG. 3: STIRAP schemes in three-level and triple-well sys- 
tems, a) Three- level scheme with pump flp{t) and Stokes 
Qsit) laser pulses and detuning A. b) Counterintuitive se- 
quence of pump and Stocks pulses with the overlap r. c) 
Triple- well scheme with the barrier couplings Qp{t) and ^ls{t) 
and detuning A. d) Three-step STIRAP transfer in the cir- 
cular well chain. The dash arrows mark three STIRAP steps. 



the adiabatic condition. 
A{t) makes ^/n'^{t) 



Toll . Indeed, the rapid rise of 



(t) large enough not only at the 
crossing point A(t) — where rt(t) is maximal but also 
at the remote flanks of il(i). As a result, the condition 
(fTO| takes place for a long time, which favors the adia- 
batic transfer. Further, as is seen from the plots d)-e), 
our scheme provides the complete transport even for the 
interacting condensate (0 < A < 2.5). The transfer is fine 
at A = 1.0 and still persists, though with fluctuations, at 
A — 2.5. For stronger interaction A — 4.0 the transfer 
breaks down. So, provided the rapid detuning change, 
our scheme ensures a robust and complete transport of 
the BEG even under a modest interaction. 

In Fig. 3 the STIRAP population transfer scheme is 
illustrated for triple-level a) and triple-well c) systems. 
Now we deal with three states (|1), |2) and |3)), pump 
and Stocks time-dependent couplings (ilp(t) and Qs{t)), 
and constant detuning A. The detuning is not crucial 
here and can be omitted 1^, 1^. Instead, the counterin- 
tuitive sequence of the couplings (Stocks precedes pump) 
and their overlap during sufficiently long time t arc im- 
portant, see plot b). Only under such conditions it is 
possible to provide a slow adiabatic evolution from the 
initial state |1) to the target state |3). The adiabatic 
condition \0{t)\ ^ 0(t) [IG,, JJi] can be written as 



nsit)np{t) - npit)nsit) 



< 1 or rir > 



(11) 



where ^Ir is the average amplitude of the pump and 
Stokes couplings during the overlap time t, d ^ ti2 — ^23 
is the relative pump-Stokes shift, and F is the width pa- 
rameter in (l7|). Following STIRAP practice, one may 
take fir ~ 0.5 AT and r w F. Then it is seen that the 
adiabatic transfer needs a strong coupling K and/or a 
long overlap time r. 
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FIG. 4; STIRAP transport of BEG in the circular well chain. 
The plots show evolution of the populations A''i (solid line), 
N2 (dash line) and (dotted line) at different values of the 
interaction and detuning, as indicated. The vertical arrows 
mark Stocks-pump coupling pairs. The coupling width and 
shift are r=3.7 and d=2. The time is dimensionless. The 
initial phase differences are zero. 



If the couplings |1) ^ |2) and |2) <-> |3) are sup- 
plemented by the additional one |3) ^ |1), then one 
may launch, via three STIRAP steps, |1) |3), |3) 
1 2), |2) ^ a cyclic adiabatic transfer, see plot 3d). 
Such a cyclic adiabatic transport can be used to generate 
the geometric Berry phase [24]. 

In Fig. 4, results of our calculations for 3 sequential 
STIRAP-like transfers (constituting altogether the cyclic 
evolution) are considered. In the plot a) we see the ro- 
bust STIRAP transfer at A = 0, i.e. without the in- 
teraction between BEC atoms and related noninearity. 
Being initially in the first weU {Ni = I, N2 = = Q), 
BEC undergoes the sequence of full population transfers, 
|1) ^ |3), |3) ^ |2), |2) ^ |1), and finally fully returns 
back to the well |1). The transfer survives at modest in- 
teraction, A < 0.4 in plots b)-c), and detuning, A < 0.2 
in plot d), but ruins at their larger values, see plot e). 
Unlike we have not found an improvement of the 
STIRAP transfer due to the detuning A. Instead, as is 
illustrated in the plot f), the detuning mainly spoils the 
transport. As compared with we consider not one 
but three STIRAP steps. Nevertheless, even in so long 
transfer chain, the transport of the interacting BEC per- 
sists at A = 0. Altogether, Fig. 4 shows that STIRAP- 
like scheme is quite reasonable for BEC transport in the 
triple- well trap 27 1. 

Note that transport in Fig. 4 is not fully adiabatic. 
Even in the best cases we see small peaks signifying a 
temporary slight populations of the intermediate states, 
hence deviation from the strict adiabatic case [3, • In 
fact we have here quasiadiabatic transfer, though with 
a high fidelity. In principle, one may improve the adi- 
abaticity upgrading the transfer parameters |24|. How- 



ever, this is not the aim of the present paper. Instead, we 
are interested in the effective and complete BEC trans- 
port between the initial and target states, regardless of 
the temporary weak population of the intermediate state. 
Moreover, our calculations show that the complete trans- 
fer can be done even at intuitive sequence of the couplings 
when the pump precedes Stokes 2J], i.e. in strictly non- 
adiabatic case. However, the adiabatic transfer has some 
advantages, e.g. it is less sensitive to the parameters of 
the process. 

In summary, we propose simple and effective adiabatic 
population transfer schemes for the complete and irre- 
versible transport of BEC in double-and triple- well traps. 
The results are straightforwardly generalized for multi- 
well traps with M > 3. The schemes work in the presence 
of a modest nonlinearity caused by interatomic interac- 
tion. Our developments open interesting perspectives for 
investigation of various geometric phases 2J|, e.g. using 
the prescription [2^. There could be also various ap- 
plications foi^Josephson qubits [l5 |. a diabatic quantum 
computing [26], topological states \\^. etc. At zero 
interaction, A = 0, our proposals are relevant for the 
transport of individual atoms. 

The work was supported by grants PVE 0067-11/2005 
(CAPES, Brazil) and 08-0200118 (RFBR, Russia). 
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